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Abstract 

New Crystal Barrel data are reported for pp — > urj and pp ojir^ir^ with to decaying to 
Tr~^7:~'K^. The wq data confirm angular distributions obtained earlier from data where u — >■ 
The new cjr/ data provide accurate measurements of vector and tensor polarisations of 
the Lo and lead to considerable improvements in masses and widths of s-channel resonances. 
A new 

JPC = 3+- J = resonance is observed with mass M = 2025 =b 20 MeV and width 
r = 145 lb 30 MeV. Polarisation is close to zero everywhere and tensor polarisations are 
large, as is the case also for pp lott^. 

Earlier publications have reported studies of pp — )■ urj [1] and pp — )■ utt^tt^ [2] for beam 
momenta 600 to 1940 MeV/c. That work concerned all-neutral final states where u decays 
to 7r°7. Those data suffer from the disadvantage that much of the information concerning u 
polarisation is lost, because it is transferred to the unmeasured polarisation of the decay photon. 
In order to recover that information, we present here new data where u — t- 7r+7r^7r°. In these 
data, the polarisation of the u is determined fully by the normal n to the decay plane of the 
u, as explained in an accompanying paper on uir^ and urjir^ final states [3]. Data for pp — )■ urj 
have also been reported by Peters [4]. 

Analyses of channels with quantum numbers C = +1, 1 = and 1 [5-7] have led to a spectrum 
of s-channel resonances consistent with almost all of the expected qq states in the mass range 
1900-2400 MeV. The data presented here lead to major improvements in the identification of 
resonances with I = 0, C = —1. All the expected states are observed except for two J^'" = 1 
states. In that sector there is a problem separating "^Si and '^Di partial waves. 

The present work follows closely the procedures described in the accompanying paper on 
I = 1, C = —1, where data for final states un^, urjir^ and 7r^7r+ are discussed. Most experimental 
details are common to that work, so we refer to it for description of the experimental set-up and 
many technical details. 

The essential experimental problem is to minimise background from r/Tr+Tr^vr^ in the corj 
data and from n^ir'n^Tr^iT^ in the un^ir^ data. An initial selection requires that both charged 
particles are produced with centre of mass angle 6, \ cos6\ < 0.65, in order to avoid edge-effects 
in the drift chamber. Here 6 is the lab angle of the u with respect to the beam. At least 11 
digitisations are required in this chamber, with at least one hit in the first three layers and at 
least one in the last three. In the preliminary selection, a 4C kinematic fit is required with 
confidence level > 5% for ?77r+7r~7r° or n'^TT^TT^TT^ir^. 

After this initial selection, clear signals are observed for u — > tt+tt^tt'^ in both sets of data, 
as shown in Figs. 1(a) and (b); all 7r° combinations are included in the latter figure. The 
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background under the u is fitted to a quadratic function of mass and the uj is fitted to a 
Gaussian. Events are selected for further processing in the mass range 760-804 MeV; the few 
events where more than one 7r"'"7r~7r'^ combination hes in this interval are rejected. Surviving 
events are then subjected to a kinematic fit to final states urj or cjtt^tt^, setting the u) mass to 
781.95 MeV. 

Table 1: Numbers of selected events (including background) 



Beam momentum 

(MeV/c) 


OUT] 


a;7r"7r" 


600 


1139 


5446 


900 


4768 


25071 


1200 


2106 


11510 


1525 


612 


4893 


1642 


1209 


11031 


1940 


469 


5241 



The final selection of ur] and wtt^tt^ requires a confidence level (CL) > 10% and greater than 
that of any background channel. As a minor refinement to check that the u is well reconstructed, 
it is required that CL{ur]) > 0.8 x CL{'K^'K^'K^ri) and > 0.5 x CL(7r+7r^47). Corresponding 
cuts are applied in the selection of cjtt^tt^ events. The efficiency with which events pass the 
final kinematic fit is determined as a function of 7r"'"7r~7r° mass, and this is used to evaluate the 
resulting background under the lo. For lot] it increases steadily with beam momentum over the 
range 13 to 16%; for c^Tr^Tr" it is 28-35%. Here we also allow for the contribution to background 
from 'wrong' combinations of the spectator tt" with tt+tt". Numbers of selected events are shown 
in Table 1. 

As explained in the accompanying paper [3], the background may be estimated in a second 
way. Decays of the ui are enhanced near the edge of its Dalitz plot because of the momentum 
dependence of the matrix element for uj decay. Figs. 1(c) and (d) show plots of the number 
of events against the square of this matrix element. One sees straight lines with intercepts 
which provide another estimate of the background; it agrees with the first within errors. This 
background is included into the partial wave analysis described below, using Monte Carlo events 
which pass the data selection; they are generated according to 7r+7r~7r°?7 or 7r+7r~7r°7r°7r° phase 
space. 

We now compare angular distributions for present ut] data with those from all-neutral data 
where u — )■ 7r*^7. Figs. l(e)-(j) show error corridors through the latter data. Points with errors 
are superposed from present data. (Both sets of data are corrected for acceptance). The absolute 
efficiency for tracking charged particles has a significant uncertainty; it is sensitive to precise 
cuts on the number of layers and the for the fit to a helix. Therefore, the absolute scale for 
present charged data is normalised to that for all-neutral data. 

Particles are thrown forwards in the lab system by the Lorentz boost due to beam momentum. 
In consequence, ui are detected efficiently only in the backward hemisphere in the centre of 
mass. This is adequate, since conservation of C-parity demands that the production angular 
distribution is symmetric forwards and backwards with respect to the beam. Figs. l(e)-(j) show 
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Figure 1: M{7i^7i~7i^) from the preliminary data selection for (a) ur] data, (b) wtt^tt" at 900 
MeV/c; the number of events v. the matrix element squared for u decay to (c) cur], (d) cutt^tt"; 
(e)-(j): curves show the error corridor for our] differential cross sections where cu decays to 7r°7; 
corrections have been applied for angular acceptance; points with errors show new results for oj 
decays to 7r"'"7r~7r°; the vertical dashed lines show the cut-off in cos^^j which has been used. 
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Figure 2: Vector polarisation Py at (a) 900, (b) 1200, (c) 1525 and (d) 1940 MeV/c compared 
with the partial wave fit (histogram) ; (e)-(h) Re T21 at the same momenta. 



that the shapes of angular distributions agree well for the backward hemisphere between the 
two sets of data. This agreement, seen particularly clearly at 900 MeV/c where statistics are 
highest, is a valuable cross-check on experimental techniques. In the forward hemisphere, we 
reject events where the acceptance for the ui drops rapidly. This requires a cut cos < 0.7 at 
600 and 900 MeV/c, cos^^^ < 0.2 at higher momenta. These cuts reject only ~ 10% of selected 
events. 

Vector and tensor polarisations of the u are determined following the procedures described 
in the accompanying paper [3]. We discuss first results for the two-body final state curj. Fig. 2 
shows values of vector polarisation Py and also Re T21 at four momenta; Fig. 3 shows values of 
Re T22 and T20. Tensor polarisations are large. Experimental values of Im T21 and Im T22 are 
consistent with zero, as predicted theoretically. 

We turn now to the partial wave analysis. This follows precisely the lines described for / = 1, 
C — —1 [1-3]. In outline, partial wave amplitudes are described as a sum 

. ^ giexp{i(j)i)BL{q)Bi{p) 

J ^ Mf-s-iM,T, ^' 

over s-channel resonances of constant width; fitted parameters are masses and widths, coupling 
constants gi and phases 0j. Blatt-Weisskopf centrifugal barrier factors Bf^ and are included for 
production with orbital angular momentum i in the pp channel and decay with orbital angular 
momentum L; p and q arc centre of mass momenta in pp and meson channels respectively. We 
adopt a radius of 0.83 fm for the centrifugal barrier radius in all partial waves up to angular 
momentum 3, as determined in Ref. [5]; this radius is increased to 1.1 fm for higher partial 
wfives. 

The spectroscopic notation is described in detail in the accompanying paper [3]. All triplet 
partial waves may couple in either the pp entrance channel or in the decay to orbital angular 
momentum L = J ± 1. A resonance is expected to have the same phase (p iii these coupled 
channels, through rescattering between the two channels. For present data there is an accurate 
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Figure 3: (a)-(d) Re T22 and (e)-(h) T20 at 900, 1200, 1525 and 1940 MeV/c, compared with 
the partial wave fit (histograms). 



determination of the relative phase only for the two 2 states; it is —5.3 ± 17.6° for the lower 
one and 1.7 ± 17.6° for the second. These are both consistent with zero, as expected. The same 
result was found also in Refs. 3 and 5, with somewhat better errors. Couphng constants for 
coupled channels are therefore fitted to a real ratio r = gL=j+i/ gL=j_i. 

Table 2 shows parameters of fitted resonances. The information from the polarisation of the 
u leads to major improvements in the ut] channel compared with the earlier work of Ref. [1]. 
For this channel, all resonances are now well determined in mass and width, except for the 
3 state at 2285 MeV. Columns 5 and 7 of Table 2 show changes in log likehhood when each 
resonance is omitted from the fit and remaining resonances are re-optimised. For convenience, 
columns 6 and 8 show corresponding values from the earlier analyses. One sees immediately a 
considerable increase in the significance of most resonances. 

Analysis of the (jjtt^tt^ channel gives less precise results for several reasons. One is that 
the uj polarisation information is distributed over 3-body phase space, though it is all used in 
the maximum likelihood fit. The decay channels which are significant are a'/2(1270), oua and 
6i(1235)7r, see Ref. [2]; there is also a contribution from fo{1500)uj at the highest two momenta. 
The bin channel is weak. The ua channel is strong for several partial waves, particuarly J^'-'^ = 
1 ; here a stands for the tth S-wave amplitude, as parametrised by Bugg, Sarantsev and Zou 
[8]. There is a problem concerning the treatment of the cua amplitude. In tttt elastic scattering, 
there is an Adler zero close to threshold, as a consequence of the nearly massless pion. In 
coupling to heavy channels, pp and a;, it is not clear whether the Adler zero should be present in 
the amplitude or not. Some examples are known where the Adler zero is definitely absent, e.g. 
J/^' — )■ ua [9]. We have explored both possibilities for every partial wave, choosing the better 
alternative for every resonance. This leads to considerable but unavoidable flexibility in the fit. 

A related problem is that there are large interferences between cua and ujf2- The precise 
phases of these amplitudes affect masses of fitted s-channel resonances; masses can shift with 
fitted phases in such a way as to keep changes in log likelihood small. It leads to rather large 
errors in resonance parameters for those partial waves where ucr contributions are big. A final 
consideration is that background is fairly high for cutt^tt^ data. Although the 3-body data 
demand large contributions to many partial waves, requiring the presence of resonances, they 
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Table 2: Resonance parameters from a combined fit to oorj and a;7r°7r°, using both uj — )■ 7r°7 
and UJ — )■ 7r+7r~7r° decays. Values in parentheses are fixed. Values of r are ratios of coupling 
constants gj+i/gj^i. Columns 5 and 7 show changes in S =log likelihood when each resonance 
is omitted from this fit and others are re-optimised; columns 6 and 8 show a comparison with 
previous work. 





Mass M 


Width r 


r 


AS 


Previous 


AS 


Previous 




(MeV) 


(MeV) 




(urj) 


AS{uj7]) 


(cutttt) 


AS{u!7l7l) 


1+- 


1965 ± 45 


345 ± 75 


-0.20 ±0.17 


1808 


209 


700 


143 


1+- 


2215 ± 40 


325 ± 55 


4.45 ± 1.16 


272 


76 


1020 


360 


3+- 


2025 ± 20 


145 ± 30 


0.60 ±0.49 


341 




234 




3+- 


2275 ± 25 


190 ±45 


-0.74 ±0.60 


253 


17 


140 


240 


1— 


1960 ± 25 


195 ± 60 


2.4 ±0.45 


1258 


281 


393 


1515 


1— 


2205 ± 30 


350 ± 90 


-3.0 ± 1.8 


180 


53 


1410 


1021 


2" 


1975 ± 20 


175 ± 25 


0.60 ±0.10 


948 


34 


1203 


478 


2" 


2195 ± 30 


225 ± 40 


0.28 ±0.59 


160 


94 


356 


590 


3— 


1945 ± 20 


115 ±22 


0.0 ± 1.0 


752 


230 


5793 


595 


3— 


2285 ± 60 


230 ± 40 


1.4 ± 1.0 


98 


46 


1156 


510 


3— 


2255 ± 15 


175 ± 30 


- 50 


145 


187 


1436 


695 


4— 


2250 ± 30 


150 ±50 


(0) 


70 






35 


5— 


~ 2250 


320 ± 95 


(0) 


677 




244 





lead to rather imprecise determinations of resonance parameters in many cases. 

We now discuss individual partial waves, beginning with singlet states. Their intensities in 
the OUT) data as a function of mass are shown in Fig. 4. Intensities of the strong contributions 
to unn data are shown in Fig. 5. There are some changes compared with Figs. 5 and 6 of Ref. 
2. In several cases, these changes arise directly from better determinations of r parameters, by 
switching of amphtudes between different L values in the final state for a given J^. 

The spectroscopic notation of Fig. 5 is illustrated by that for panel (b) as an example. There, 
= '^^~^^£j denotes the initial pp state, where i — 2 combines with spin s = 1 to make J — 2; 
P = (-1)^ and C = (-1)^+*. The J = 2 states may decay with L = 0, 2 or 4. In almost all 
cases, only the lowest L is significant for cutt^tt" data, because of the centrifugal barrier. In Fig. 
5(b), ^5*2 = ^'^ '^^Lj denotes the final state ujf2 with s' = 2, L = 0, J = 2. It is usually necessary 
to include all available s' values for the lowest L. Examples are Figs. 5(h) and (i), where s' = 1 
and 2 respectively, combining with L = 1 to make J — 1; Clebsch- Gordon coefficients tend to 
favour larger values of s'. 

For J-^*^ = 3"'"", the present data reveal a new resonance which escaped detection in the 
earlier analyses. It lies at 2025 ± 20 MeV with P = 145 ± 30 MeV. A qq state is expected 
close to this mass as the singlet partner to ^^2, ^F^ and states observed clearly in Ref. [5]. 
Spin-spin splitting is believed to come from a S function at zero radius [10]. In the absence of 
such sphtting in F-waves, the singlet state is expected to he at the centroid of the triplet states 
weighted by their multiplicities, namely at 2024 ± 5 MeV; this agrees with the present result. 

A second 3"*"" state is observed at 2275 ± 25 MeV. This compares with the prediction of 
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Figure 4: Intensities of partial waves fitted to cor] data; L are orbital angular momenta in the 
decay, and 3S1, for example, denotes pp ^Si. 



2292 ± 9 MeV from triplet states. Both of these states are determined very largely by the 
OUT) data. Although they make sizeable contributions to a;7r°7r° data, errors are large there for 
masses and widths. 

For J^*^ = 1^* , two states arc observed at 1965 and 2215 McV. The accuracy of the mass 
determination of the lower state comes mostly from urj, although there is a large contribution in 
un^n^ too. The error is sizeable, because this resonance is at the bottom of the available mass 
range. The upper state is required by the strong ^Pi(a;/2) contribution of Fig. 5(i). Both masses 
are slightly lower than those from Rcfs. [1] and [2], as a consequence of the new u polarisation 
information. The qq centrifugal barrier is weaker for 1+ states than for 3^, so it is to be expected 
that 1+ will resonate somewhat lower than S"*", as is observed. The lower l+~ state now appears 
in ujTj almost purely with L — decays; that is a large change from the earlier result of Ref . [1] . 
The new polarisation data for the uj have provided a major improvement in the determination of 
parameters r = gi^j_^_i/gi^^j_i, which describe the ratio of coupling constants g for L = J ± 1. 

We turn now to triplet partial waves, beginning with the 1 sector, which is the most 
difficult. For J-^*-^ = 1 , there are large contributions in cutt^tt" and these data are mostly 
responsible for fixing the upper resonance at 2205 MeV. The our) data do however produce the 
best determination of the lowest 1 resonance at 1960 MeV; it is weaker in un^n^, though 
clearly visible in the 6i(1235)7r channel. Fig. 5(a). It couples strongly to ^Di, where its mass is 
best determined; it is probably the radial excitation of w(1650). 

Four 1 states are expected in this mass range: two ^Si and two ^Di. Unfortunately, ^5*1 
and ^Di are not well separated for remaining 1 states in the absence of data from a polarised 
target. From a comparison with 7 = 1, C = — 1 [3] where data are available from a polarised 
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Figure 5: Intensities of the large partial waves fitted to cutt^tt" data. The spectroscopic notation 
is discussed in the text. 
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target for the final state tt'tt"^, higher resonances are expected around 2110-2150 and 2230 
MeV and at a higher mass 2350-2400 MeV. It is possible that the state we report at 2205 
MeV is a bhirred combination of two states in the mass range 2110-2230 MeV. It has a large 
uj(T contribution, shown in Figs. 5(m) and (n), and this contribution is somewhat flexible, as 
discussed above. This resonance does have quite a large ^Di amphtude, and is consistent with 
the state expected around 2230 MeV. We have attempted to put two 1 resonances into 
the mass range 2100-2230 MeV, but they collapse to a single state. 

In earlier work of Refs. [1] and [2], there was evidence for a further 1 state at 2300 ± 45 
MeV. Although an improvement in log likelihood of ~ 300 is possible by including this state, 
there is no well defined optimum for its mass and width. Therefore it is omitted from the final 
fit. Omitting it has little effect on other partial waves. 

For = 2", two resonances are definitely required for two reasons. Firstly, the phase 
variation visible in the Argand diagram of Fig. 7 below for L = 1 is larger than can be 
provided by a single resonance. Secondly, the structure observed in Fig. 5(b) requires two 
resonances. The lower one is well determined by both ur] data (M = 1981 ± 23 MeV) and wtttt 
data (1973 ± 24 MeV). The upper one is determined better by urf data. By comparison with 
/ = 1,C = — lof Ref. 3, it is expected around 2235 MeV. In present data it appears somewhat 
lower, at 2195 ± 30 MeV, though within the combined errors. 

For J^'^ = 3 , three states are definitely required. The lowest contributes a huge f20J signal 
in 7r°7r°a; data, see Fig. 5(d); it is also quite large in urj. This state is expected to lie close to 
the very well defined I — 1 state p3(1982). That resonance is defined by extensive pp — )■ 7r~7r+ 
data at 100 MeV/c steps of momentum down to 360 McV/c (a mass of 1910 MeV) [11]; its 
mass is 1982 ± 14 MeV and the width 188 ± 24 MeV. In present data, the / = state optimises 
at 1944 ± 16 MeV in cun^n^ and at 1951 ± 21 MeV in ur]. Table 2 quotes an average value 
1945 ± 20 MeV; the error allows for the small discrepancy between channels. It is remarkably 
narrow in both data sets: F = 115 ± 22 MeV. This narrow width is required to fit a very rapid 
variation in differential cross sections between beam momenta of 600 and 900 MeV/c. However, 
we warn that this state lies right at the bottom of the available mass range, so there could be 
some systematic error in determining its mass and width. 

A ^Gz state is required strongly by both ur] and un'^n^ data. It optimises at M = 2255 ± 15 
MeV with F = 175 ± 30 MeV. It is conspicuous because of its different helicity components from 
^Da, hence different tensor polarisations. A ^D^ state is also required at a similar mass. However, 
this is the least well determined of all the resonances. It makes a fairly small contribution to uf] 
data. For cutt^tt^, a strong contribution is required. However, the mass drifts up to 2310 MeV 
with only a small improvement in log likelihood; in the earlier analysis of Ref. 2, it optimised 
instead at 2180 ±40 MeV. These instabilities in mass are correlated with flexibility in fltting the 
oja component in cutttt. For the present solution, our flnal best estimate for the mass (mostly 
from cur] data) is 2285 ± 60 MeV. For both ur] and cutt^tt^, this state is obscured by the very 
large low mass contributions from C(;3(1945). 

A small but well determined ^G^ peak is now required in our]; a ^G^ state is required strongly 
by both Lur] and a;7r°7r° data. For all G-states, the intensity peaks strongly at ~ 2300 MeV, 
see Figs. 4 and 5. However, there is a very strong centrifugal barrier in the pp channel. In 
consequence, we flnd anomalously low resonance masses for all three G-states, compared with 
neighbouring F-states, which cluster from 2260 to 2305 MeV. We use standard Blatt-Weisskopf 
centrifugal barrier factors, which are derived assuming a square well potential. This may not be 
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Figure 6: A comparison of for resonances with straight-line trajectories against radial exci- 
tation number n; the slope of 1.143 GeV^ is taken from Ref. [5]. In (b), the ^D^ trajectory is 
displaced one place to the right in n in order to resolve it from ^Di. 

a good approximation for very high partial waves such as G-states. 

Fig. 6 shows a plot of observed resonances versus M^, where M is mass. They conform 
closely to straight-line trajectories resembling Regge trajectories, to which they are related. 
They are compared in the figure with a slope of 1.143 GeV^. That slope is determined from 
very precise results for 7 = C = -|-1 [5]; for those quantum numbers, data are available from 
seven channels including valuable 7r~7r"'" data from a polarised target. Within the errors of 
present determinations, there is good agreement with this slope. 

Fig. 7 shows Argand diagrams for all partial waves. A remarkable feature of the vector 
polarisations of Fig. 2 for ut] is that they are everywhere small. The same feature has been 
observed for a;7r°. Vector polarisation depends on the imaginary parts of interferences between 
partial waves. The small observed polarisation requires coherence amongst all partial waves. It 
is discussed in detail in the accompanying paper. 

In summary, the new data provide a considerable improvement in parameters of several 
resonances. This arises from more accurate polarisation information for the u, leading to better 
phase determinations. A new but expected 3+~ resonance is observed at 2025 ± 20 MeV with 
r = 145 ± 30 MeV. 
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Figure 7: Argand diagrams for fitted ujt] partial waves. Crosses show beam momenta 600, 900, 
1050, 1200, 1350, 1525, 1642, 1800 and 1940 MeV/c; all move anti-clockwise with increasing 
beam momentum. 
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